B -> h + h'- modes at CDF 
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We review CDF results and prospects on decays of B hadrons in two charged charmless hadrons. 
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I. INTRODUCTION 

The rich production of all types of B hadrons at the 
Tevatron Collider at Fermilab allows studying charmless 
two-body decays both in known and new modes. Using 
a sample of J Cdt ~ 1 fb _1 of data, CDF performed a 
search for all possible modes of neutral bottom hadrons 
in two charged charmless hadrons (p, K or w) 

The CDFII detector is a multipurpose magnetic spec- 
trometer surrounded by calorimeters and muon detec- 
tors 0. A silicon micro-strip detector (SVXII) and a 
cylindrical drift chamber (COT) immersed in a 1.4 T 
solenoidal magnetic field reconstruct charged particles 
within pseudorapidity \r)\ < 1.0, with a transverse mo- 
mentum resolution <j pt /pt — 0.15%p*r/(GeV/c). This 
yields a mass resolution of ~ 22 MeV/c 2 for B — > h + h'~ 
decays, which is an important ingredient of this analysis. 
Particle identification information (PID) is obtained from 
the specific energy loss by ionization (dE/dx) of charged 
particles in the COT. This yields a nearly-constant sep- 
aration of 1.4 a between pions and kaons of momenta 
>2 GeV/c. A three-level trigger sytem allows selecting 
events by requiring the presence of at least two tracks 
with large impact parameters relative to the beam axis. 



II. DATA SAMPLE AND RECONSTRUCTION 

B hadron candidates are initially selected by forming 
pairs of oppositely-charged tracks with > 2 GeV/c, 
transverse opening-angle 20° < A0 < 135°, and pr(l) + 
Pt(2) > 5.5 GeV/c. In addition, both tracks are required 
to have a large transverse impact parameter do relative to 
the pp interaction vertex (100 /im < d a < 1 mm). The B 
candidate is required to point back to the primary vertex 
(do(B) < 140 /im), and to have travelled a transverse 
distance L xy {B) > 200 /xm. 

Most of the above cuts are implemented at the trigger 
level. In the offline analysis, additional cuts are imposed 
on isolation Ib and the quality of the fit (x 2 ) to the 
3D decay vertex of the B hadron candidate. 

Final selection cuts are determined by an optimization 
procedure, based on minimizing the expected uncertainty 
of the physics observables to be measured. Two differ- 
ent sets of cuts are used, optimized respectively for best 



resolution on Arp(B° — ► K + ir~) (loose cuts), and for 
best sensitivity [J] for the discovery of the yet unobserved 
B® —> K~tt + mode (tight cuts). The looser set of cuts 
is also used for measuring the decay rates of the largest- 
yields modes, while the tighter set is used for the other 
rare modes. 

The invariant mass distribution of the candidates, with 
an arbitrary pion mass assignment to both tracks, shows 
a single large peak in the B mass range, formed by several 
overlapping modes (Fig. [JJ. 



III. MEASUREMENT METHODOLOGY 

The different modes are statistically separated and in- 
dividually measured by means of an unbinned maximum- 
Likelihood fit, combining kinematics and PID. Kinematic 
information is summarized by three loosely correlated ob- 
servables: (a) the mass M W7r calculated with the pion 
mass assignment to both particles; (b) the signed mo- 
mentum imbalance a = (1 — Pi/p2)q\, where p\ fj^) is 
the lower (higher) of the particle momenta, and q\ is the 
sign of the charge of the particle of momentum p\\ (c) 
the scalar sum of particle momenta p to t — Pi + P2- The 
above variables allow evaluating the mass of the B candi- 
date for any mass assigment to the decay products. PID 
information is given by a dE/dx measurement for each 
track. The Likelihood for the i th event is then: 
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where index 'A(E)' labels the physics (combinatorial) 
background-related quantities, and index j runs over 
the twelve possible B — > h + h'~ and A° — > ph~ modes 
and conjugates having distinguishable final states (e.g. 

B° -> K+tj- and B -► K~ 7r + are distinct, while 
B° — > 7r + 7r - and B —> ir + ir~ are treated as a single 
component). The fj are the signal fractions to be deter- 
mined by the fit, together with the background fraction 
parameters b and /a- 

The shape of the mass distribution of each single chan- 
nel accounts for non-Gaussian tails, both from resolution 
and from emission of photons in the final state, which is 
simulated on the basis of analytical QED calculations 0] ■ 
The quality of this model was checked on a large sample 
(~500k) of D° —> K~it + decays. The mass distribution 
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FIG. 1: Invariant mass distribution of B? s \ — » h + h candidates passing the loose (left) or tight (right) selection cuts. The 
pion mass is assigned to both tracks. Cumulative projections of the likelihood fit for each mode are overlaid. 



of the combinatorial background is fit to a smooth func- 
tion, while the physics background is parameterized by 
an 'Argus function' Q smeared with our mass resolu- 
tion. Kinematical distributions for the signal are repre- 
sented by analytical expressions, while for the combina- 
torial background are parameterized from the mass side- 
bands of data • 

The dE/dx response of the detector to kaons and pions 
was calibrated from a sample of 1.5M D* + — > D°ir + — > 
[K~ir + ]ir + decays, where the D° decay products are 
identified by the charge of the D* + pion. The PID term 
for the background allows for pion, kaon, proton, and 
electron components, which are free to vary indepen- 
dently in the fit. Background muons are indistinguish- 
able from pions with the available dE/dx resolution and 
are therefore included as a part of the pion component. 

To avoid the large uncertainties associated to produc- 
tion cross sections and reconstruction efficiency, branch- 
ing fractions are measured relative to the B° — > K + ir~ 
mode, and then normalized to the world-average value 
of B(B° -> K+n-) 8]. Upper limits [| are quoted for 
modes in which no significant signal is observed. 

To convert the raw signal fractions returned by the fit 
into relative branching fractions, corrections are applied 
for different efficiencies of trigger and offline selection re- 
quirements for different decay modes. Corrections re- 
lated to decay kinematics are determined from the detec- 
tor simulation, while others are measured on data using 
control samples. The dominant contributions to the sys- 
tematic uncertainty come from: statistical uncertainty on 
isolation efficiency ratio (for B® modes); uncertainty on 
the dE/dx calibration and parameterization; and uncer- 
tainty on the combinatorial background model. Smaller 
systematic uncertainties are assigned for: trigger efficien- 
cies; physics background shape and kinematics; B meson 
masses and lifetimes. 



IV. RESULTS 

A. Rare Modes 

The search for rare modes is performed using the 'tight' 
selection. The fit allows for the presence of any com- 
ponent of the form B — > h + h'~ or A° — > ph~ where 
h, h' = K or 7r, with the yield as a free parameter. Final 
results are reported in Table Q] where fa and f s indicate 
the production fractions respectively of B° and B® from 
fragmentation of a b quark in pp collisions. 

The results provide the first observation of the B® — > 
K~ir + mode, with a significance of 8.2a, which in- 
cludes systematic uncertainties and is evaluated from 
Monte Carlo samples of background without signal. The 
branching fraction of this mode is significantly sensi- 
tive to the value of angle 7; however, predictions ob- 
tained from different methods differ. Our measurement 
B(B° K-TT+) = (5.0 ± 0.75 ± 1.0) x 10~ 6 is in agree- 
ment with the prediction in [Toj , but is lower than most 
other predictions [lj, El, EH . 

No evidence is found for modes B® — > 7r + 7r~ or 
B° — > K + K~, in agreement with expectations of signif- 
icantly smaller branching fractions. The measurement 
B(B° -> K+R-) = (0.39 ± 0.16 ± 0.12) x 10~ 6 has 
the same precision of the other current measurements 
from T(4S I ) experiments 0], although the upper limit is 
weaker due to the observed positive central value. The 
B® — > 7r + 7r~ upper limit is improved with respect to the 
previous best limit, obtained from the analysis of a sub- 
sample of the present data [l4[ . The sensitivity to both 
B° — » K + K~ and B® — » 7r + 7r~ is now close to the upper 
end of theoretically expected range [III 0) EH (3 , and 
it will be interesting to see the results from the larger 
samples being accumulated. These modes proceed only 
through annihilation and exchange diagrams, which are 
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currently poorly known and a source of significant un- 
certainty in many theoretical calculations. Measuring 
of constraining both these channels is particularly use- 
ful since the physics parameters governing the strength 
of peng uin-annihilation can be extracted from their ra- 
tio [rzf- 

In the same sample, we also get to observe charmless 
decays of a B baryon for the first time: A" — > pir~ (6a) 
and A" — » pK~ (11.5<r). We measure the ratio of branch- 
ing fractions of these modes as B(A® — > pn~)/B(A° — > 
pK~) = 0.66 ± 0.14 ± 0.08, in good agreement with the 
expected range [0.60, 0.62] from [HI]. Work is in progress 
towards a measurement of individual branching fractions 
and CP asymmmetries, that are expected to be non- 
negligible and are sensitive to possible SUSY contribu- 
tions fl9j . 



B. CP asymmetries 

We can measure from our data the CP asymmetries 
of both B° and B® decays in the self-tagging final state 
JiT ± 7r T . The asymmetry of the B® mode is measured 
with the tight selection used in previous section, while 
the looser selection is used for the B° mode. 

The raw asymmetries returned by the fitting procedure 
need to be corrected for possible detector and procedural 
biases. This has been done by measuring the asymmetry 
in a sample of 1M prompt D° — » K~tt + decays, recon- 
structed and selected with the same computer code used 
for analysis of the B — > h + h'~ sample and similar selec- 
tion criteria. Given the smallness of the CP asymmetry 
expected in the D° — > K~ir + mode (<< 1%), the exper- 
imentally measured asymmetry provides a good determi- 
nation of the measurement bias, including asymmetries 
in the dE/dx response or other possible unanticipated ef- 
fects. The observed effect (0.6 ± 0.14)% is in good agree- 
ment with the expected K + / K~ asymmetry due to the 
different probability of interaction with the detector ma- 
terial. 

The result A CP {B° -> K+tt-) = -0.086±0.023±0.009 
is in agreement with the world-average [H, and is the 
second most precise measurement. The updated world 
average A^ e (B -> K+tt-) = -0.095 ± 0.013 has a 
significance of 7a (previously 6a). Comparison of this 
average with the asymmetry in the similar mode of 
the B + , shows a deviation of 4.8<r (previously 4.6c). 
While it has been argued in the past that the asymme- 
tries in these two modes should be equal in the stan- 
dard model due to isospin symmetry, this is not any- 
more considered a reliable test [2l|, [23, [23| . Conversely, 
it has been argued 0, HH that a much more robust 
test of the Standard Model origin of the asymmetry 
of the B° — > K+ir~ mode is the comparison with the 
corresponding asymmetry in the B® —* K~it+ mode, 
where the final state is identical, thus canceling pos- 
sible effects of final state interactions. The predicted 
equality of rate differences T(B — > K~ir+) — T(B° —> 



K+tt-) = T{B° S -> AT-7T+) - r(B° s -> K+n-) is very 
robust under the Standard Model, while it would be 
completely fortuitous under a New Physics scenario, be- 
cause it is produced by interference of very different am- 
plitudes in the two cases. In addition, the smallness 
of the B(B° S -> K-TT+) (Tab. fl}, makes the predicted 
asymmetry large (~ 40%) and therefore experimentally 
more accessible. Using our tight set of cuts, we find 
A C p(B° s -> K-TT+) = 0.39 ± 0.15 ± 0.08. This value fa- 
vors the large CP asymmetry predicted by the Standard 
Model and has the correct sign, but is still compatible 
with zero (significance just above 2a). By combining our 
measurement with the world-average for the B° we ob- 

tain - -i — ^=n : - = 0.84 ± 0.42 ± 0.15, m 

r(B°->K-n+)-r(Ws-> K+tt-) ' 

agreement with the Standard Model expectation of 1.0. 
It will be very interesting to see if this agreement per- 
sists with more data. Given the large expected asymme- 
try, the SM asymmetry will be visible as a 3<r effect with 
1.5 fb _1 of data, and a 5a effect before reaching 4 fb _1 . 
A non observation of this asymmetry would indicate a 
non-SM source of CP violation. 



C. Precision Branching Fractions 

The sample selection for Aqp(B q — > K+tt~) also pro- 
vides good measurements of the branching fractions of 
the 'large yield' modes B® — > K + K~ and B° — * 7r + 7r~ '. 
We obtain B(B° S -> K+K~) = (24.4 ± 1.4 ± 4,6) x 10" 6 , 
in agreement with current predictions [jH l23l [25l [26| and 
with the previous CDF measurement [lj], although the 
lower central value now indicates a smaller U-spin break- 
ing effect. Work is in progress to reduce the systematic 
uncertainty which dominates the resolution of the present 
preliminary results. 

A substantial yield is also available in the B° — > 7r + 7r~ 
mode, allowing a measurement of the branching fraction: 
B(B° -> 7r+7r-)= (5.10 ± 0.33 ± 0.36) x 10~ 6 . This has 
the same precision, and is in agreement with the current 
results from T(45) [|. 



V. PROSPECTS 

CDF results based on 1 fb _1 are beginning to show the 
Tevatron potential for B physics. The sample on tape is 
now almost doubled, and the luminosity keeps increas- 
ing; the default plan is to collect 8 fb _1 by year 2009. 
Most of the systematic uncertainty in the measurements 
described above is determined by statistical uncertain- 
ties in the calibration samples, and is therefore expected 
that the precision will keep increasing with statistics. 
Highlights from the full runll samples will be: Aqp in 
B° — > K + tt~ at 1% level; 5-sigma observation of direct 
Aqp in B® — » K~ir + (or alternatively the discovery of 
non-SM CP violation) ; first Aqp measurements in the 
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TABLE I: Results from the loose (top) and tight (bottom) selections. Absolute branching fractions are normalized to the the 
world-average values B(B° -> K + tv~) = (19.7±0.6) x 10~ 6 and / s = (10.4±1.4)% and f d = (39.8±1.0)% 0. The first quoted 
uncertainty is statistical, the second is systematic. Limits are at 90% CL. 
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1121 ± 63 

1307 ± 64 
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B° - 


-> K+K 




■* pK~ 


A2- 
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230 ± 34 ± 16 

26 ± 16 ± 14 
61 ± 25 ± 35 
156 ± 20 ± 11 
110 ± 18 ± 16 



Quantity 



Measurement 



A CP 

B(S"-»K+7r- 
7J" B(B°->K+7r-) 



-0.086 ± 0.023 ± 0.009 
0.259 ± 0.017 ± 0.016 
0.324 ± 0.019 ± 0.041 



S B(S»- 



A CP 

f s B(B°^: 
fd B(B°-,K+-k-) 
B(b"~>K+ K~) 



■ K+tt-) 



B{B®~- 

B(A°- 



K+ir- 
■fm~) 



B(A»^pK-) 



0.066 ± 0.010 ± 0.010 

0.39 ± 0.15 ± 0.08 
0.007 ± 0.004 ± 0.005 
0.020 ± 0.008 ± 0.006 
0.66 ± 0.14 ± 0.08 



e(io- b ) 



5.10 ± 0.33 ± 0.36 
24.4 ± 1.4 ± 4.6 



5.0 ± 0.75 ± 1.0 

0.53 ± 0.31 ± 0.40 (< 1.36) 
0.39 ± 0.16 ± 0.12 (< 0.7) 



A" ; and tight constraint, or even observation, of annihi- 
lation modes. 

In addition to the above, time-dependent measure- 
ments will be performed. Resolutions on time-dependent 
CP asymmetries can be predicted from the yields, flavor 
tagger perfomance (e£> 2 ), and effective S/B. Tagger per- 
formance was optimized to eD 2 — 5.3% for the B® mixing 
analysis; we assume a performance eD 2 = 4% can be ob- 
tained for the B° with similar methods. For the other 
parameters we assume the values of current analysis with 
no further improvements. 

The CP asymmetries in the B° — > 7r + 7r~ mode can be 
measured with a resolution <jacp = 0.15 with 6 fb _1 , 



which will offer an interesting additional measurement of 
similar precision to the currently available results from 
T(45) , that still show a disagreement. 

The performance in the measurement of CP asymmetry 
in the B® — > K + K~ mode depends on the assumed effec- 
tive S/B (~ 1, depending on the selection), and proper- 
time resolution (between 70 and 100 fs -1 ). A resolution 
in the range [0.15, 0.3] is expected for both direct and 
mixing Aqp from a sample of 6 fb _1 . These asymmetries 
are related in the SM to the B° asymmetries and angle 
7 via U-spin [23], [2?J , and are sensitive to possible SUSY 
effects [26]. In addition, it is possible to measure the B° s 
lifetime in a CP eigenstate, see 28] for details. 



C-conjugate modes are implied and branching fractions 
intended as CP-averages, unless otherwise stated. 
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